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Abstract 



We calculate the nucleon parameters in nuclear matter using the QCD sum rules 
approach in gas approximation. Terms up to in the operator product expan- 
sion (OPE) are taken into account. The higher moments of the nucleon structure 
functions are included. The complete set of the nucleon expectation values of the 
four-quark operators is employed. Earlier the lack of information on these values 
has been the main obstacle for the further development of the approach. We show 
that the values of the four-quark condensates are consistent with the assumptions 
about the convergence of the OPE. Inclusion of these condensates and of the non- 
locality of the vector condensate are important for the calculation of the nucleon 
parameters. The nucleon vector self-energy S^, and the nucleon effective mass m* 
are expressed in terms of the in-medium values of QCD condensates. The numerical 
results for these parameters at the saturation value of the density agree with those 
obtained by the methods of nuclear physics. 



1 Introduction 



The QCD sum rules were invented by Shifman et al. PP to express the hadron parameters 
through the vacuum expectation values of QCD operators. Being initially used for the 
mesons, the method was expanded by loffe fS] to the description of the baryons. The 
approach succeeded in describing the static characteristics as well as some of the dynamical 
characteristics of the hadrons in vacuum — see, e.g., the reviews jHlll!- 
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The basic idea is to consider the correlation function no(g^) describing the propagation 
of the system with the quantum numbers of the hadron, in the different regions of values of 
the momentum g^, where certain information on its behavior is available. The asymptotic 
freedom of QCD enables to present Ilo{q^) at — > — oo as a power series in and the 
QCD coupling as- On the other hand, the imaginary part of Ilo{q^) at > can be 
described in terms of the observable hadrons. This prompts to consider the dispersion 
relation for the function no(g^) PP 



at g^ — * — oo. The coefficients of the expansion of the left-hand side (Ihs) of the function 
no(g^) in powers of g~^ are the expectation values of the local operators constructed of the 
quark and gluon fields, which are called " condensates" . Such presentation, known as the 
operator product expansion (OPE) [S] provides the perturbative expansion of the short- 
distance effects, while the nonperturbative physics is contained in the condensates. The 
usual treatment of the right-hand side (rhs) of Eq. (1) consists in "pole -|- continuum" 
presentation, in which the lowest lying pole is singled out while the higher states are 
approximated by the continuum. Thus, Eq. (1) ties the values of QCD condensates with 
the characteristics of the lowest hadronic state. Such interpretation requires that the 
contribution of the pole to the rhs of Eq. (1) exceeds the contribution of the continuum. 

The OPE of the Ihs of Eq. (1) becomes increasingly valid, when the value of |g^| 
increases. On the other hand, the "pole + continuum" model becomes more accurate 
when |g^| decreases. The important assumption is that the two presentations are close in 
a certain intermediate region of the values of g^. To improve the overlap of the QCD and 
the phenomenological descriptions, one usually applies certain mathematical tools, i.e. 
the Borel transform. The Borel transformed dispersion relations (1) are known as QCD 
sum rules [H |2] . 

For example, the QCD sum rules for the nucleon provided a connection between the 
nucleon mass and the scalar quark condensate (0|gg|0) [T. Similar relations have been 
obtained for the magnetic moments of the nucleons [6 , etc. 

Later the QCD sum rules were applied for the investigation of modified nucleon pa- 
rameters in nuclear matter [71 |H] . They were based on the Borel-transformed dispersion 
relation for the function Ilm{q) describing the propagation of the system with the quan- 
tum numbers of the nucleon (the proton) in nuclear matter. Considering nuclear matter 
as a system of A nucleons with momenta pi, one introduces the vector 



which is thus p ~ (m, 0) in the rest frame of the matter. The function Ilm{q) can 
be presented as Ilm{q) = ^m{q'^,'^{p,q)) with the arbitrary function (p{p,q) being kept 
constant in the dispersion relations in g^. 




(1) 
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The spectrum of the function Ilm{(l) is much more comphcated than that of the func- 
tion no(g^). The choice of the function (p{p,q) is dictated by attempts to separate the 
singularities connected with the nucleon in the matter from those connected with the 
properties of the matter itself. Since the latter manifest themselves as singularities in the 
variable s = {p + g)^, the separation can be done by putting (p{p,q) = {p + qY and by 
fixing 13 IHl H 

V?(p,g) = {p + qf = s = AEIj, (3) 
with Eqf being the relativistic value of the nucleon energy at the Fermi surface. 
The general form of the function Urn can thus be presented as 

n™(g) = q,l^K.{ci\ s) + p,7^^^(g^ s) + ml{q\ s) . (4) 

The in-medium QCD sum rules are the Borel-transformed dispersion relations for the 
components Iii^{q^,s) (j = q,p,I) 

It was shown in [3 |Hl El CD! that the "pole + continuum" model for the rhs of Eq. (5) can 
be used at least until we do not include the higher order terms of the density expansions 
of the functions Ili^{q'^,s). Thus, one can expect that the characteristics of the nucleon 
in nuclear matter can be expressed through the in-medium values of QCD condensate. 

In the lowest order of OPE the problem was approached in [3 |S1 El HO] • It was noticed 
that the condensates of the lowest dimensions {d = 3, 4) can either be calculated or ex- 
pressed through the observables. The vector condensate Vfj,{p) = (M| Z^i ^*(0)7m?'(0)I-^) 
is proportional to the density of the matter p, being 

v^ip) = VN,P , VN, = {N\ ^ g^(0)7^g^(0)|iV) . (6) 

i 

Here the upper index i denotes the quark flavor. In the rest frame of the matter we get 
Vfi{p) = v{p)6^o, VNf, = v^Sf^o with 

v{p) = vnP ; vn = S (6') 

being just the number of the valence quarks in the nucleon. The scalar condensate is 

kM = {M\Y^q\0)q\Q)\M) = Ko + ^ip) , 

i 

Ko = kUO), k{p) = KmP + ---, KN = {N\J2q\0)q\0)\N) . (7) 

i 

Here the dots denote the terms which are nonlinear in p. The expectation value (A^| J2i (fq^\N) 
is related to the nN sigma term ct^at, i.e. |TT] 

= {N\uu + dd\N) = Jfll^ (8) 
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with rriu^d standing for the current masses of the hght quarks. 

Turning to the condensates of dimension d = 4, we find for the gluon condensate [7j 

gM = {M\^G\0)\M) = go + 9{p), 

71 

9o = 9miO) , gip) = 9nP + ■■■ (9) 
with the nucleon expectation value 

= {Nl^G'mN) ^ --m, (10) 
71 y 

obtained in [T^ in a model-independent way. 

Also, the nonlocal condensate (A^|g(0)7^g(x)|A^) provides the contributions of d = 4 
and those of the higher dimension. The term of the dimension d = 4 is 

{N\q^{Oh,D,q\0)\N) = ^g^^-^^mx, (11) 

with X2 standing for the second moment of the nucleon structure function [U] . The nonlo- 
cality of the scalar operator q{Q)q{x) manifests itself in the higher terms of the operator 
expansion. The nonlocality of the product of the gluon operators is not expected to be 
important because of the minor contribution of the gluon expectation value to the nucleon 
parameters. 

The shift of the position of the nucleon pole, which in the linear approximation can 
be identified with the single-particle potential energy of the nucleon, was expressed as a 
linear combination of the condensates of the lowest dimension |7, .8] . The vector and scalar 
expectation values appeared to be the most important ingredients. Their contributions 
cancelled to large extent, reproducing the familiar features of the Walecka model [TSj . 
An alternative approach was developed in |^ with the dispersion relations in the 
time component go at three-dimensional momentum |q| being fixed. It provided a similar 
result. 

The lack of knowledge about the in-medium expectation values of the higher dimension 
became the obstacle for the development of both approaches. One of such expectation 
values is the scalar four-quark condensate {M\qqqq\M) . It was noticed in jHJ HH] that 
the configuration 2{Q\qq\Gi) {N\qq\N) p (with p standing for the baryon density) is one of 
those, which composed the in-medium expectation value of the operator qqqq. In the gas 
approximation the expectation value of the colour-singlet operator is 

{M\qqqq\M) = (0|gggg|0) + 2p (0|gg|0) (Ar|gg|Ar) + p {N\{qqqq)int\N) (12) 

with the last term describing the " internal" action of the operators inside the nucleon. In 
the "ground-state saturation approximation" (also called "factorization approximation") 
formulated in the last term of the rhs of Eq. (12) vanishes. This would lead to the 
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change {M\qqqq\M) — {0\qqqq\0) = 2p(0|gg|0)(A^|gg|A^) of the value of the scalar four- 
quark condensate. Assuming this approximation one would be forced to conclude that 
the four-quark scalar condensate plays the crucial role in QCD sum rules, causing doubts 
on the convergence of OPE. The numerical results would contradict the known nuclear 
phenomenology [16|. 

There have been some attempts to get rid of the contribution of the four-quark con- 
densates, applying the differential operators [H] or by choosing the form of the in-medium 
function Ilm{<l) which does not couple to the four-quark condensates |T7j. However, some 
of the information appeared to be lost in the former case, while there still remained some 
unknown condensates in the latter case. Anyway, there was no consistent analysis of sum 
rules with the inclusion of the four-quark condensates until now. On the other hand, there 
are some indications that the second term of the rhs of Eq. (12) does not provide the true 
scale for the in-medium modification of the value of the scalar four-quark condensates. 
The calculations, carried out in |^ predicted strong cancellation between the second and 
the third terms in the rhs of Eq. (12). Also, the arguments based on chiral counting and 
supporting the violation of the in- medium factorization have been given in [T^ . 

In the present paper we build and solve the QCD sum rules in nuclear matter in 
the gas approximation with the account of the condensates up to the dimension d = 6. 
This means that we include the terms of the order 1/g^ of the OPE (recall that the 
leading OPE terms are of the order g^lng^). This requires the inclusion of the four-quark 
condensates uT^uuT'^u, dT^ ddT"^ d and uT^udT'^ d with F"^'^ standing for the basic 
4x4 matrices, corresponding to the scalar, pseudoscalar, vector, pseudovector (axial) 
and tensor structures. 

In the gas approximation the in-medium expectation value of any operator A is 

{M\A\M) = (0|i|0) +p(iV|i|A^) (13) 

with |A^) standing for the state vector of the free unpolarized nucleon. Since we include 
only terms linear in p, we can neglect the Fermi motion of the nucleons of the matter. 
Thus we put 

s = 4m^ (14) 

in Eq. (3). Having in mind the future extension of the approach, we shall keep the 
dependence on s, using Eq. (14) for the specific computations. 

We consider symmetric nuclear matter with an equal density of the protons and neu- 
trons pp = pn = p/2. 

The nucleon expectation values of the lowest dimensions can either be calculated in 
a model-independent way or expressed through the observables. The calculations of the 
four-quark condensates require model assumptions on the structure of the nucleon. The 
complete set of the four-quark condensates was obtained in [201 by using features of the 
perturbative chiral quark model (PCQM). The chiral quark model, originally suggested 
in [21], was developed further in |22|- In the PCQM the nucleon is treated as a system 
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of relativistic valence quarks moving in an effective static field. The valence quarks are 
supplemented by a perturbative cloud of pseudoscalar mesons, in agreement with the 
requirements of the chiral symmetry. In pUj a simple version close to the SU(2) flavor 
PCQM, which includes only the pions, has been used. 

There are three types of contributions to the four-quark condensate in the framework 
of this approach. All four operators can act on the valence quarks. Also, four operators 
can act on the pion. There is also a possibility that two of the operators act on the valence 
quarks while the other two act on the pions. Following [20] we speak of the "interference 
terms" in the latter case. 

To obtain the contribution of the pion cloud, we need the expectation values of the 
four-quark operators in pions. The latter have been deduced in j2Sl by using the current 
algebra technique. We obtain a remarkable cancellation of the pion contributions in the 
function Ilm{q)- This cancellation takes place in any model of the nucleon which treats the 
pion cloud perturbatively. Thus, the contributions of the four-quark condensates come 
from the terms, determined by the valence quarks only and from the interference terms. 

We find the contribution of the four-quark condensate to the in-medium modification 
Um — Ho of the function IIo to be much smaller, then one could expect by assuming 
the "in-medium factorization approximation". These terms are about 4-5 times smaller 
than the leading ones of the OPE series. This is consistent with the hypothesis of the 
convergence of OPE. 

Thus, we obtain three sum rules equations for the functions n^(g^, s), n^(g^, s) and 
n^(g^,s) introduced in Eq.(4). The in- medium characteristics of the nucleon, i.e. the 
vector self-energy S^, and the effective mass m* are the unknowns of these equations. 
There are two unknown parameters more, i.e. the residue at the nucleon pole X*^ and the 
continuum threshold W^. All these characteristics will be obtained from the QCD sum 
rules. 

The dependence of the rhs of Eq. (5) on the parameters, which are expected to 
be determined, is not linear (except the dependence on Aj^). Thus, even in the gas 
approximation the behavior of these parameters with p is linear only at sufficiently small 
values of the density. We consider two approaches to the problem. In the linearized 
case we determine only the linear parts of the in-medium modifications of the hadron 
(nucleon) and continuum parameters. We construct the combination of the sum rules for 
the function 11^ — Hq in such a way, that two of the equations determine the values of 
S^, and m* — m separately. Note that it is possible to write the equation in which the 
parameter m*—m is the only unknown, only because the proton has a definite space parity. 
The third equation enables to find the in- medium changes of the parameters 6X^ = Aj^ — Ag 
and 6W'^ = — Wq. We express the nucleon characteristics and m* — m in terms of 
the vector, scalar, gluon and four-quark condensates and of the moments of the structure 
functions. 

In the nonlinearized version we do not assume the in-medium changes of the parame- 
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ters to be small. The three equations for the Borel transformed function 11^(5') — Ho^q"^) 
enable to obtain the values of T,y, \*^ and W^. At densities p of the order of the saturation 
value po the values of S^, and m* — m appear to coincide within 25% and 10% accuracy 
with the values provided by the linear version. This causes somewhat larger difference in 
the values of the potential energy U{p) which still has reasonable values. 

Inclusion of the four-quark condensates and of the higher moments of the structure 
functions diminish the OPE value of the nucleon vector self-energy S^, by about 25% each. 
As to the scalar self-energy m* — m, the four-quark condensates provide contribution of 
the same order as the leading OPE term. However, this contribution is almost totally 
compensated by the account of the higher moments of the structure functions. Thus, the 
value of the m* — m is very close to that given by the leading order of OPE. 

The structure of the paper is as follows. In Sec. II we present the sum rules in a form, 
which is convenient for our analysis. In Sec. Ill and IV we calculate the contribution 
of the dimension d = 6, i.e. the expansion of the nucleon structure functions and the 
four-quark condensates. In Sec. V - VII we present the solutions in the linearized and 
nonlinearized forms. We discuss and summarize the results in Sec. VIII and IX. 

2 General equations 
2.1 Sum rules in vacuum 

To make the paper self-consistent, we recall the main points of the QCD sum rules ap- 
proach in vacuum mi2]- The function UqIq'^) (often referred to as "polarization operator") 
is presented as 



with j being the three-quark local operator (often referred to as "current") with the proton 
quantum numbers. The usual choice is 



where T denotes a transpose and C is the charge conjugation matrix. The upper indices 
denote the colors. 

The Ihs of Eq. (1) is approximated by several lowest terms of OPE, i.e. Ilo{q^) ~ 
n^^^(g^). The empirical data are used for the spectral function Im Ilo{q^) on the rhs 
of Eq. (1). Namely, it is known, that the lowest lying state is the bound state of three 
quarks, which manifests itself as a pole in the (unknown) point A;^ = m^. Since the next 
singularity is the branching point k"^ = Wpf^ = (m + m^)^, one can present 




(15) 



(16) 



Imno(A;2) = \l5{k^ 



m')+fikyie-W') 



(17) 
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with being the residue at the nucleon pole. Thus, Eq. (1) takes the form 

^rV) = ;^ + i/^*^ (18) 

pn 

Of course, the detailed structure of the spectral density /(A;^) cannot be resolved in such 
an approach. The further approximations are based on the asymptotic behavior 

m = ^ Anr^(A:^) (19) 

at fc^ 3> \q'^\ with A denoting the discontinuity. The discontinuity is caused by the 
logarithmic contributions of the perturbative OPE terms. The usual ansatz consists in 
extrapolation of Eq. (19) to all the values of k"^, replacing also the physical threshold W^^^ 
by the unknown effective threshold Wq, i.e. 



1 7 m , _ 1 7 ^^o'^^ik 



2^ 



TT J k^ — q"^ 27ri J k"^ — q 

pn 

Thus Eq. (1) takes the form 



dk'' ^ I "V" ' . (20) 



tOPE/„2\ ^"n , /" An^^^(^^) ,,2 



no"""(?^) - + 7^ 7 dk' ■ (21) 

— q^ 2m J — q^ 

WS 

The Ihs of Eq. (21) contains QCD condensates. The rhs of Eq. (21) contains three 
unknown parameters: m, and Wq. The OPE becomes increasingly true when the value 
\q'^\ increases. The "pole + continuum" model is more accurate at the smaller values of 
|g^|. Thus one can expect Eq. (21) to be true in a certain limited interval of the values 
of |g^|. To improve the overlap of the OPE and the phenomenological description one 
usually applies the Borel transform defined as 



g2 = -g2. = Q^n (22) 



with M called Borel mass. It is important in the applications to the sum rules that the 
Borel transform eliminates the polynomials and emphasizes the contribution of the lowest 
state in rhs of Eq. (21) due to the relation 



B -^—^ = e-'/^' . (23) 
— q^ 

The Borel-transformed form of Eq. (21) reads 

n^^^(M2) = A^e— '/^' + -^ / dA;2e-^'/^'An^^^(fc2) (24) 

27rz J 
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and is known as QCD sum rules. Actually, there are two sum rules for the structures Uq 
and IIq of the function no(g) = g^7^no(q'^) +IIlQ{q^) with / standing for the unit matrix. 

It appeared to be more convenient to work with Eq. (24) multiplied by the numerical 
factor 32-7?^. The two sum rules for the nucleon in vacuum can be presented in the form 

m 

Ll{M\W^) = Ao(M2) . (25) 
mM^,W^) = mAo(M2) (26) 



with 



Here = 32n^X%, 



Ao(M2) = A^e-'"'/*'' . (27) 



Lg(M2, 1^2)^ 32^4 /n^,,OPi^^^2)_^ f rfA;V=VM2Ang'^^^(A;' 
V 27ri J 

CO 

Ll{M\W^) = 327r^fno^'^^^(M2)-— / dk^-''" /^lil''''''' 



:28) 



The Ihs of Eqs. (25) and (26) jSE] have been obtained by including of the condensates 
of dimension ci = 8, i.e. with the account of the terms of the order l/q^ in OPE of the 
functions IIo '^'^^ 

r , n r,, . ah 111 a, o 1 , , 

Ll{M\Wl) = ^^M'E,-- + — -^a'— (30) 

with the traditional notations a = — (27r)^(0|gg|0) = —2tx'^hq (we assumed the isotopic 
invariance (0|um|0) = (0|Jd|0) = (0|gg|0)), h = (2TifgQ, /xg = 0.8 GeV^. Here Ei are the 
functions of the ratio W^/M^: Ei = EiiWi/M'^). They are given by the formulas 

Eo(a;) = 1 - e-", Ei{x) = I - {I + x)e-\ E2(x) = 1 - + x + y j e"". (31) 

The factor 

,,_2x InMVA^ ^ ^ 

accounts for the anomalous dimension, i.e. the most important corrections of the or- 
der as enhanced by the "large logarithms". In Eq. (32) A = Aqcd = 0.15 GeV, 
while u = 0.5 GeV is the normalization point of the characteristic involved. Note that 
the two last terms on the rhs of Eq. (29) originate from the four-quark condensates 
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{Q\uT^uuT^u\0) and can be expressed through the single term ((0|gg|0))^ only in frame- 
work of the factorization hypothesis mi2]- Also, the last term on the rhs of Eq. (30) is 
the six-quark condensate, evaluated in the same approximation. 

The matching of the Ihs and rhs of Eqs. (25) and (26) have been achieved [21 El in the 
domain 

0.8 GeV^ < < 1.4 GeV^ (33) 
providing the values of the vacuum parameters 

= 1.9 GeV^ ; = 2.2 GeV^ (33') 

if m = 0.94 GeV. 



2.2 Sum rules in nuclear matter 

The OPE terms of the polarization operator in nuclear matter 

Um{q) = I I d^xe'^^-^ {M\Tj{x)jiO)\M) (34) 



contains the in-medium values of QCD condensates. Some of these condensates vanish in 
the vacuum, obtaining non-zero values only in the medium. The other ones just change 
their values compared to the vacuum ones. 

The spectrum of the function Ilm.{(l) is much more complicated, than that of the 
vacuum function no(g^). However, [H |H1 13 UD] the spectrum of the function Ilm{q'^, s) 
at fixed value of s can be described by the "pole + continuum" model at least until we 
include the terms of the order in the OPE of Umiq"^, s). 

The description of the nucleon pole is based on the general expression for the propa- 
gator 

Gn' = {G%y'-^ (35) 
with G% = [q^'y^ — m)~^ being the propagator of the free nucleon, while 

E = g^7'^S, + - p^7^Sp + (36) 
m 

is the general form of the self-energy of the nucleon in nuclear mater. In the kinematics, 
determined by Eq. (3) we obtain 

^ ^ qf,!'' -Pf,l''{K/m) +m* 

q^ -mi, 

with 
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The new position of the nucleon pole is 



2 (s — m^)Sj,/m — + m*^ 



while 



< = ^ \ly , " , (39) 

1 + L^/m 



(40) 



:i-S,)(l + S,/m) 



Thus, we shall present the dispersion relations for the functions n^(g^, s) {i = q,p, I) 
determined by Eq. (4) in the form 



with bg = 1, bp = —Hy, bj = m*. The Borel-transformed sum rules take the form 

LliMMV^) = AUM'); (42) 
Ll{M\Wl) = -E^AUM'); (43) 
Li{M\W^) = m*AUM') (44) 

with 

Am{M') = A^^e-"^™/^' . (45) 

Here 

X*^ = Xl-Z (46) 
is the effective value of the residue in nuclear matter. 
We present the Ihs of Eqs. as 

with tjn{M'^,W^) standing for the lowest order OPE terms, u\^{M'^,W^) denoting the 
contribution of the higher moments of the structure functions, while uj''^{M'^) provides the 
contribution of the four-quark condensates. We write, correspondingly, Lq = £q + uOq for 
the Ihs of the vacuum sum rules presented by Eqs. (25) and (26). We present also 

t{M\wl) = e^{M\wl) - e,{M\wl) , 

uj\M') = u;UM') - uji^iM') . (48) 
In these notations the lowest order OPE terms are 

^0 = ^4/9 + 4 14/9 ' ^0 = 0, io = 2aM Ei^ (49) 
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and 



= aM\Wl)v{p) + fl{M\Wl)g{p) , 
P = nM\Wl)v{p), 

(50) 

with 

Jv 3 ^^4/9 ' h ^4/9 

Iv 3 ^4/9 

/^^ = -47r2M^Ei„ . (51) 

The functions f (p), ^(p), (^(p) are determined by Eqs. (6'), (7), (9). The notation E^m 
{k = 0, 1, 2) means that the functions depend on the ratio W^/M^. Actually, the higher 
moments of the structure functions of the nucleon have been neglected in Eqs. ()5U|) and 

(EH). 



3 Accounting for x-dependence of the operators. 
Contributions of the higher moments and 
of the higher twists of the structure functions 

The calculation of the function Ilm{q'^, s) defined by Eq. (34) is based on the presentation 
of the single-quark propagator in the medium 

(M|Tgj,(x)g^(0)|M) = G,^(a:)-l(M|g^(0)7^g^(a:)|M)7^^- J(M|r(0)g*(x)|M)5„^ (52) 

with G{x) = (zx^7^)/(27r^x'^) being the free propagator of the quark in the chiral limit. 
Recall that i denotes the light quark flavor. In the lowest orders of OPE two of the 
quarks are described by the free propagators and only one of the quarks is presented by 
the second or the third term of the rhs of Eq. ()52|) . 

At X = the matrix elements in the second and third terms on the rhs are just the 
vector and scalar condensates defined by Eqs. (6) and (7). The contribution of the bilocal 
configurations can be expressed in terms of the higher moments and twists of the nucleon 
structure functions 0. 

The bilocal operators on the rhs of Eq. ()52|) are not gauge invariant. The gauge 
invariant expression, achieved by substitution |24j 

q'{x) = q\0)+XaD^q'{0) + ^XaXfsD^Dpq'{0) + --- (53) 
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with Da standing for the covariant derivatives, provides the infinite set of the local con- 
densates. The expectation values depend on the variables (px) and x^. In the gas approx- 
imation we only need the nucleon matrix elements Eq. (13). For the vector structure the 
general form is 

e;ix) = {N\q\0)j,q\x)\N) = ^ <Pli{px) , x') + tx .m^Hpx) , x') (54) 

with g*(a;) defined by Eq. ()53|). 

Expansion in powers of corresponds to the expansion of the function Ilm{(l) in 
powers of g^. To obtain the terms of the order it is sufficient to include two lowest 
terms of the expansions in powers of x^. One can present 0123 

1 

</>l(,)((pa;),x2) = /ciae-^'^(^-)/:(,)(a,x2) (55) 



with ^ 

/a(6)(«,a;2) = ^a(6)(a) + ga;^"^^C(fe)(«)- (56) 

Here ?7a(a) = /„(«, 0) is the contribution of the quarks with the fiavor i to the asymptotics 
of the nucleon structure function 77(a) = ?7a(a) + (a), normalized by the condition 

1 

dari{a) = 3 (57) 



with the rhs presenting just the number of the valence quarks in the nucleon. Thus, 
expansion of the function ^p\{px) = 0Jj((px), 0) in powers of (px) is expressed through the 
moments of the distributions r)1{c()- The moments are well known — at least, those, which 
are numerically important. Also, the first moment of the distribution ^a(tt) = C('^)+^a('^) 



J. 

e = / ^ -0.3 (5^ 



was calculated in 120] by QCD sum rules method. The moments of the function rjKa) can 
be obtained by using the equations of motion Da'j°'q^{x) = miq^{x). Thus, in the chiral 
limit P 

1 /, . o, 1 



(^i") = 5(R«')-4(e)) 



6 

Here we denoted 



(a) = • (59) 



(/)= I'dafia) (60) 
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for any function /(a). 

Note that the nonlocahty of the scalar condensate, i.e. of the last term on the rhs of 
Eq. (|52p does not manifest itself in the terms up to l/q"^. The first derivative in (px), as 
well as all the derivatives of the odd order vanish in the chiral limit due to QCD equation 
of motion. The next to leading order of the expansion in powers of x"^ vanishes due to 
certain cancellations [5| as well as in the case of vacuum pi for the particular choice of 
the operator j{x) presented by Eq. (16). We do not account for the nonlocality of the 
gluon operators, since the gluon expectation values play the minor role in our sum rules. 

Now we are ready to calculate the contributions n„;(g) of the nonlocal vector conden- 
sate to the polarization operator Ilm{q). We express n„; in terms of the proton expectation 
values 6'^(x) = {p\(t{0)'~f^q'^{x)\p) . Employing the isotopic invariance we obtain 

^ni{q) = -J^( x'^-^ + x{x, r + e')] e^(''^) • p, (61) 



contributing to the vector structures q and p of the polarization operator limiq)- Here we 
denoted 

a = a^i^. (62) 

Using Eq.(j3D) we obtain 

Ilni{q)=Ki{<D+Ki{<l) (63) 

with 

= 5 / ^ (^^P^^ + S:{xpm + e^^"^^ ■ p (64) 

n!;.(g) = -5/^x(0;; + 0^)e^(-).p 

We present each of the terms H^f ■* as the sum + H^"^'''* , corresponding to the two 

terms of the expansion in powers of x"^ in Eq. lj^Uj) . In particular, the contribution 11^°, 
which is numerically most important, can be presented as: 



.'2 ^ I '2 



6m7r^ 



J daq'{pq')ln{-^)r]a{a) + J daq'^ ln{-^)r]a{a) 

^ ^ 



P (65) 



with rjaia) = Va{(^) + vt{(^)^ q' = q — pa (see Appendix A). From Eqs. (3) and (14) one 
finds (pq) = (s — — q^)/2. The cutoff Ac will be eliminated by the Borel transform. 

Presenting g'^ = —{l + a){Q'^ + A^) where Q"^ = — g^, A'^{a) = a{s — m'^ — m'^a)/{l + a) 
we see that the second term of the expression Ing'^ = Ing^ + In(^) does not have a cut, 
running to infinity, but has a finite cut. This singularity requires a special treatment 
in QCD sum rules. On the other hand, it is the singularity in the u channel of the 
interaction of the baryon current with the quark of the nucleon of matter. It corresponds 
to the exchange terms on the rhs of the sum rules. In this paper we neglect the nonlocal 
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singularities, thus claiming for the description of the nucleon in the Hartree approximation. 
However, we include the regular smooth dependence on the higher moments. 

The contributions 11^^ and H^i^''^ can be expressed in terms of the moments of the 
functions rjl and ^* ^ (see Appendix A). Since the higher moments of the functions ?7*(tt), 
as well as the value of ^ are small, we include only the lowest moments of the functions 
ril{a) and the first moments of the functions C(q^) — Eq. (j^ . The last of the equalities 
enables us to neglect the contribution of the functions ^^(a). 

Finally, the higher moments and higher twists of the nucleon structure functions pro- 
vide the contributions to the Ihs of the sum rules — Eq. 
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,2 . 5 



in- 



(66) 



in- 



3L4/9 



-5(M^Ei„ - (s - m^)M'Eom){vc^) 



Here we denote L = L{M'^). Parameter ^ is defined by Eq. (I58|) . 



4 Contribution of the four-quark condensates 

The four-quark expectation values contribute to the OPE terms 1/g^ of the function 
Hm(5')- Now only one quark is determined by the free propagator Gq{x). Two other 
quarks are described by the last term of the two-quark propagator 

{M\Tq^{x)qp{Q)qp{x)qM\M) = [G,{x)f - ^-{M\qV'' q\M)G,{x)V^p 
- i(M|gr''g|M)G,(a;)rJ, + ^mq^"" dd^"^ (l\M)T^pT^^, (67) 

with F'''"'^ being the basic 4x4 matrices 

F^ = /, F^^ = 75, F^ = 7;., r^ = 7M75, = ^(7^7. - 7.7^), (68) 

acting on the Lorentz indices of the quark operators. Equation (j67j) is analogous to Eq. 
(I52|) for the single-quark propagator. We did not display the color indices in Eq. (jHTj), 
keeping in mind that the quark operators are color antisymmetric — Eq. (16). One can 
write an equation similar to Eq. ()67|) for the quarks of different flavors. 
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Introducing the notations 



iy^^(p) = {M\uT''uuT^u\M); R^^ {p) = {M\dT''duT^u\M) (69) 
we write in the gas approximation 

H^^'^ip) = H^^'^iO) + p/.^^; R^'^ip) = i?^^(0) + pr^^ . (70) 
The characteristics h^"^ and r^^ can be presented as 

/i^^ = ^[{0\uT''u\0){N\uT^u\N) + {0\uT^u\0){N\uT''u\N)) 

+ {N\{uT''u-uT^u)ir.t\N) ; (71) 

r^^ = ^[{0\dT''d\0){N\uT^u\N) + {0\uT^u\0){N\dT''d\N)) 

+ {N\{dT''d-uT''u)int\N) . (72) 

Here the lower index "int" means that all the four operators are acting inside the nucleon. 
The coefficients 5/6 and 2/3 on rhs of Eqs. (f7T|) and (f?^ present the weights of the color- 
antisymmetric states — see Appendix B. These equations are consistent with Eq. (13) if 
we assume that some of the single-particle operators which compose the operator A can 
act on the vacuum state vector — see also [TUj . 

The contribution of the four-quark expectation values to the in-medium change of the 
polarization operator can be written as 



(n)4, = (n„)4, - (no)4, = ^ E Pxrh'''' + E r^yr^^ • (73) 

^ \X,Y X,Y J 

Here pxY and txy are certain matrices in Dirac space. They can be obtained by using 
the general expression for the function Ilm{(l) presented in ^Hj 

PXY = ^ Tr (7ar^7^r^)757"g7^75 ; 
16 

TXY = ^ Tr (7„g7^r^)757"rV75 , q = q.^ ■ (74) 

Here = 1 if has a vector or tensor structure, while = — 1 in the scalar, pseu- 
doscalar and axial cases. The sign is determined by that of the commutator between 
matrix and the charge conjugation matrix C — Eq. (16). 

The products pxYh^^ obtain nonzero values if the matrices F^ and F^ have the 
same Lorentz structure. In this case all the structures presented by Eq. (jH^ contribute 
to {Jlm)iq- The products txy^^^ do not turn to zero only if F^ has a vector or axial 
structure. In the latter case F"^ should be an axial matrix as well. In the former case F"^ 
can be either Lorentz scalar or Lorentz vector. 
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We denote h^^ = , fixx = f^x and r^^ = , txx = ^x for the similar Lorentz 
structures X and Y. The scalar and pseudoscalar expectation values are Lorentz scalars. 
Thus, their contributions can be expressed through single parameters. The latter is true 
also for the scalar- vector expectation value r^^ . We obtain 

= -| , f^Ps = ^ , = (75) 

In the other channels the four-quark condensates have more complicated structure. In 
the vector and axial channels 



UV(A) _ V{A) .V{A)m^ 

T = ar^^^g,. + br^^^Pl^.. (76) 



Using Eqs. (HU, ((721) we obtain 



^yh =-a^q-bh jJAh =af^q + bj, — (77) 



and 



In the tensor channel 
with 



rvr^ = (-10ar-26r)g-26r^, 

= (-6af-2bf)q + 2bf^. (7J 



tpu,pT = (ppPp9uT + PuPrgtip - PpPrQup " PuPpQpr ] (80) 



and 



The complete set of the four-quark expectation values ,b^ , b^ was obtained in 
PU] by using the approach motivated by the perturbative chiral quark model (PCQM) 
[7Tll22j. As explained in Introduction, the valence quarks are treated as the relativistic 
constituent quarks, while the sea quarks are approximated by those of perturbatively 
treated pions. 

There are three types of contributions to the expectation values in the approach of 
PU] . All four operators can act on the constituent quarks. Also, four operators can act 
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on the pions. There are also the "interference terms" with two of the operators acting on 
the valence quarks while the other two act on the pions. 

The contribution, corresponding to all four operators acting on pions is expressed in 
terms of the pion expectation values of the four-quark operators. The distribution of the 
pion field is determined by the PCQM. The contribution is 

{TiAq)p^ans = \Y.{^l''\^lxuT'' uuV'' u + ArxdV"" duV'' u\t:'^) \ ^ (82) 
16g2 j dml 

with S standing for the sum of the self-energy and pion-exchange contributions, while "a" 
denotes the pion isotopic states. Using the values of the four-quark operators averaged 
over pions [2Sj, we find that 

^/ix(vr"|Mr^MMr^M|vr") +4^r_x(vr"|(Jr^rfMr^M|7r") = . (83) 

Due to Eq. (jH!^ we can omit the contributions to the second terms of the rhs of 
Eqs. fl7T|) and (f?^ which are caused by the pions only. Since the terms (0|gg|0) (7r|gg|7r) 
emerge as the ingredients of the expectation values (vr|gggg|7r) [21], the cancellation (jHHj) 
influences the first terms of rhs of Eqs. (frT|) and (f7^ as well. Thus, in order to calculate 
the rhs of Eq. (j73j) it is sufficient to substitute for the operators with the same flavor 

h"" = 2- -{<d\uT''u\{)){N\{uT''uWN) + {N\{uT''uuT''u)i\N) . (84) 
6 

Here the lower index "f" means that the operators act on the valence quarks only. The 
lower index "1" corresponds to the sum of the term in which all the four operators act on 
the valence quarks and the term in which two of the operators act on the valence quarks 
while the other two act on pions. Of course, the first term in the rhs of Eq. ()84|) obtains 
a nonvanishing value only in the scalar case T-^ = I. 

The expectation values of the operators of different fiavors, providing nonvanishing 
contributions to the rhs of Eq. ()72j) are the scalar- vector condensate 

= 2 ■ '^{0\dd\0){N\u^^u\N) + {N\{ddu^^u)i\N) (85) 

and 

X 



r 



{N\{dT^duT^u)^\N) (86) 

with X standing for vector or axial structures. In the first term in the rhs of Eq. (j85p the 
nonlocality of the vector condensate is included. 

The meaning of the lower index "1" is the same as in Eq. ()84p . 

Using the complete set of the nucleon four-quark expectation values [20], we obtain 

(n)4. = [aiI + + Ai^^r^ (87) 
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with the coefficients 

A\ = 0.25 Alg = -0.57 Ai^ = 1.90 (88) 
and with the conventional notation 

a = -(27r)^(0|MM|0). (89) 

We use the value (0|mm|0) = (-241 MeV)^, corresponding to a=0.55 GeV^, employed in 0. 
Note that a is just a convenient scale for presentation of the results. It does not reflect 
the chiral properties of 114^. 

We can trace the structure of the three terms, composing Il4q determined by Eq.(|H7j)- 
see Appendix C. The q term results mainly as the sum of the expectation value of the 
product of the four tt-quark operators, described by the ffist (factorized) term on the rhs 
of Eq. (j84j) ■ and that of the product of two u and two d-quark operators in the axial 
channel-Eq. fl86|) . The p term is determined mostly by the expectation value (jHBjl in the 
vector channel. The contribution proportional to the unit matrix / is determined by the 
scalar- vector expectation value (jH3j) . It is dominated by the ffist (factorized) term on the 
rhs, while the second term diminished the value by about 30%. 

The contributions of the four-quark condensates to the Ihs of the Borel transformed 
sum rules 012I)-(jll|) are 

= u'^p ; u'^ = Aijig ; 

2 

fl^ = -Svr^a, = -8n^^—a; = -Sn^ma. (90) 

Note that we can modify our model approach by employing a more sophisticated 
model for the pion. Namely, among the interference terms contributing to the four-quark 
condensates, there is so-called "vertex interference", in which one of the vertices of the 
self-energy of the valence quark is replaced by the four-quark operator. Some of such terms 
contain the matrix elements (0|u75(i|7r~) and {0\d'y5u\7i^) , contributing to the expectation 
values (A^|'U75(i(i75u| A^), being connected with the matrix elements {N\dr^ duT^ u\N) of 
all structures by the Fierz transform. On the other hand, they depend on the values of 
the quark masses, since (01^75^17?") = — ^^^^J' with M^, {F^,) denoting the mass (decay 
constant) of pion. In a somewhat straightforward approach one substitutes the current 
quark masses. Following more sophisticated models of the pions ,27j one should substitute 
the constituent quark masses, thus obtaining much smaller values. In the latter approach 

Al = - 0.11, (91) 

while the values of A^^ and A^^ remain unchanged. In this case we find the larger cancel- 
lation between the ffist term of rhs of Eq. ()84|) and the contribution coming from rhs of 
Eq. (jSB|) . The latter is dominated by the vector expectation values. 
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5 Sum rules in nuclear matter 



Actually, we shall solve the sum rules for the difference of the operators in nuclear matter 
and in vacuum: 

L\M\ Wl, Wl) = AUM') - Ao(M2) ; (92) 
L^iM^W^) = -S,A„(M=^); (93) 
L'{M\ Wl, W^) = m*A„(M2) - mK^iW) (94) 

with L*(M2, Wl, W^) = Ll{M\ Wl) - Li{M\ W^). The ingredients of Eqs. 
are defined by Eqs. (25), (26), and 

Note that we took into account the anomalous dimensions only for the leading OPE 
terms g^lng^ and Ing^, neglecting the anomalous dimensions of the OPE terms. 

Although the anomalous dimensions of the four-quark condensates are known 28J, the 
anomalous dimension matrix is not diagonal in the basis determined by Eq. ()68p. The 
calculation of this matrix in our basis is a separate work which will be presented in further 
publications. We use the nucleon structure functions presented in [21], which include the 
anomalous dimensions of the structure functions. 

We solve Eqs. ()92|l - (jM|l in the same interval of the values as it has been done in 
vacuum Eq. (33). 

Since Eqs. (I92p - (j94p are not linear, the behavior of the in-medium parameters is not 
linear in p even if we limit ourselves to the gas approximation. However, if the density p 
is small enough, we can try the linear approximation, assuming the linear dependence of 
the nucleon characteristics on the density of matter. 

6 Sum rules in the linear approximation 

It is instructive to express the density in units of the observable saturation density 

Po = 0.17fm-3 = 1.3- lO^^GeV^ . (95) 
The parameters which will be determined from the sum rules can be presented as 

= ttjj — ; m =m + as — ; OA = A„ - Ag = — ; 
Po Po Po 

SW^ = Wl-W^ = aw — . (96) 

Po 

To obtain the parameters in the linear approximation we put Z = 1 and find 

S„ = Ep , m* = m(l + Sq) + , A;^ = A^ , 

= m,(l + Sq) + S„ + = m* + S„ (97) 
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in Eqs. (38)-(40) and We put s = Am^ in Eq. (39). 

Expansion of the Ihs of Eqs. provides the equations 

(/,^(M^ W^)v^ + fl{M\ W',)gN + u%{M') + w%) po 

- I flA- (fls + aj^^ I e -aw , (98) 

(/P(M, W^)VM + muUM') + muj^j) p, = -a.A^e"'"'/*^' ; (99) 

(/i(M^ W,')K^ + c.^(M^)) PO = [a,m - (a. + a.)^) e^^/*^^ 

+asAoe -aw • UUUj 

Note that in this form all three equations are tied. One can build up another set of 
equations with the functions L^, — mL'^ and L'^ as the Ihs. In this case the unknowns 
a„ and are determined from the separate equations. The third equation determines the 
values of ax and aw- We introduce 

n{M',W^) = p,fi{M\W^)—^{k = v,g,K), 

^0 



2\ 



T^(M^) = po/L (101) 



^0 

with the functions fl and fl^ defined by Eqs. (j3T|) and (plj) . We present 

T^iM^ W^)VN + mTP^{M\ W^) + mT^Al^ = - a, (102) 

+T^{M')Ai^-mT^{M')Al = a,; (103) 
T,«(M^ W^)VN + T^^(M^ W^)gN + T„^(M^ W^) + n{M')A% 

X 2m 1 1 dmU^.W^) , 

with Aq being defined by Eq. (27). The characteristics a^ and a^ are found from Eqs. 
fll02p and ()103p and are substituted into Eq. ()104p . The latter determines the values of 
a\ and aw- 

Note that there is one more approximation in the rhs of Eq. ()102p . Namely, we 
neglected the value 
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since there is about 80% cancellation between the two terms on the Ihs in the interval 
determined by Eq.(33). This is due to the positive parity of the nucleon state — See 
Appendix D. 

The values of the QCD parameters which enter the Ihs of Eq. (|lU2j) - (jl(J4p are deter- 
mined by Eqs. (6'), (8), (10), (65) and (86). The expectation value = {N\uu + dd\N) 
is connected to the pion- nucleon sigma term a-j^N by Eq.(8). The value of a-j^N can be 
extracted from the data on low-energy ttN scattering, being expressed through the ob- 
servable S term (S^^r) [201 • The value 

a^N = (45 ± 7) MeV (105) 

corresponds to T,-„iy=64: MeV [HI]. We shall present the specific values, corresponding 
to kat = 8. This value corresponds to ct^at = 45 MeV and the sum of the light quark 
masses m„ -|- = 11 MeV. There is an uncertainty in the value of kn due to the errors 
in determination of the values of a.j,N and m„ -|- ma. We also present the dependence of 
the characteristics of the nucleon on the value of k^. 

The values of the parameters 

= 0.108 GeV , a, = -0.178 GeV , ax = -1.29 GeV^ aw = -0.81 GeV^ (106) 

are obtained by minimization of the relative difference of the rhs and Ihs of Eqs. (|102p - 
(|104p by the chi-square method. The solution of these equations is illustrated in Fig. 
1. Note, that if we construct the equation which is the difference of Eqs. ()104|) . p05p 
the function of in the Ihs should be approximated by the constant value + a^, 
having the meaning of the potential energy. Such approximation holds with much better 
accuracy then the separate Eqs. p04p and ()105|) for the self-energies. The solution p06|) 
corresponds to the values 

= 108 MeV , m* -m = -178 MeV , -j- = -0.67 , —j- = -0.37 . (107) 

Although the sets of Eqs. (j^ - pOO|) and ()102|) - ()104p are mathematically identical, a pro- 
cedure of matching of the two sides of the equations may lead to somewhat different 
solutions. Applying the same procedure of minimization to the set of Equations (j^Hjl - 
(Unni) we find a^ = 0.108 GeV a, = -0.254 GeV, ax = -1.61 GeV^ aw = -0.91 GeV^. 
Thus the parameters and SW"^ are determined with somewhat larger uncertainties 
than the self-energy Ej^. 

As we noted at the end of Sec. IV, our model approach to the calculation of the 
four-quark condensates can be modified by using more sophisticated models of the pions 
P7] . i.e. by the account of the constituent quark masses. Using the value of Al^ given by 
Eq. dnH), we obtain from Eqs. (ITn^ - (fTni|l 

= 108 MeV , m* -m = -203 MeV , —j- = -0.71 , — ^ = -0.41 (108) 

Aq W/q 
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at the saturation density p = po. Thus, this change of the value uj% results in the change 
of the nucleon parameters by less than 15%. 

Note that the functions Tj(M^) defined by Eq. (|l(Jip (j = v, g, K,u,u!;i = q,p,I) 
depend on rather weakly. Thus, approximating 

t;(m2) = q , (109) 

we can replace in the Ihs of Eqs. (jl(J2j) - (|l(J4p the functions Tj(M^) by the constant 
coefficients Cj. Numerically the most important functions T^(M^) and T^(M^) can be 
approximated by the constant values with the errors of about 4% and 8%. The largest 
errors of about 25% emerge in the averaging of the functions T^. This solves the problem 
of expressing the in-medium change of nucleon parameters through the values of the 
condensates 

= -(C>;v + mC^ + mC2A^,)^ ; (110) 



Po 



rrf-m= (cin^ - mCy^ - mClg^ - mCl + CIA{^ - mClAl^ ^ . (Ill) 

Po 

The coefficients in the rhs of Eqs. (|110|) and (jlllj) are 

C9 = -0.062, = 0.011 GeY'\ = -0.067, = -0.074, 
d = -0.042 GeV , Ci = -0.064 GeV , 

CP = -0.090 GeV , = -0.095 , = 0.094 . (112) 

Equations (jllO|) and (jlll|) reproduce the values of J2y and m* provided by Eqs. (jl07j) 
with the accuracy of 15% and 6% correspondingly. 



7 Beyond the linear approximation 

Now we do not assume the linear dependence of the nucleon parameters on the density p. 
We find the values S^, m*, Aj^ and which minimize the difference between Ihs and 
rhs of Eqs. (j92|) - (j94|) . The consistency of the Ihs and rhs is illustrated by Fig. 2. At the 
saturation density p = po we obtain 

S^ = 150MeV, m*-m = -200 MeV , A;;^ = 1.25 GeV^ = 2.11 GeVl (113) 

The two last numbers correspond to the relative shifts 6X'^/\l = —0.35 and 6W'^/Wq = 
—0.03. Thus the linear approximation is true at p ~ po with the accuracy of about 25% 
for the vector self-energy and about 10% for the scalar one. The linear approximation 
overestimates the shift of the effective threshold. 
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Recall, that we presented the numerical results for = 8. The dependence on the 
value of Kjy is shown in Fig. 3. The density dependence of the nucleon parameters at 
kn = 8 is shown in Fig. 4. 

Using Eq. (|9ip for the value of Af^ we obtain the results which are close to those 
presented by Eq. ()113|) 

E^ = 142MeV, m*-m = -223 MeV , X*^ = 1.2AGeV'^, _ 2.09 GeVl (114) 

Note that the difference between the linear and nonlinear solutions has a strong effect 
on the value of the nucleon potential energy 

f/(p) = K{p) + m*{p)-m , (115) 

which is about -40 MeV and -70 MeV for the solutions ()113p and ()114|) at the phenomeno- 
logical saturation point p = po . 

8 Discussion 

It is instructive to follow how the values of the nucleon self-energies change, while we 
include the various contributions of the Ihs of the sum rules. The solutions of the general 
equations (92)- (94) are illustrated by Fig. 5. At the saturation density po the vector 
self-energy S^, and the effective mass m* are 340 MeV and 750 MeV correspondingly, if 
only the terms P -Eq. ()50|) are included in L* of Eq. ()47p . One can see from Fig. 5 that 
the higher moments of the structure functions and the four-quark condensates subtract 
about 100 Mev each from the value of S^. On the contrary, the two contributions to m* 
cancel to large extent, with the four-quark condensate subtracting about 200 MeV, and 
the moments of the structure functions adding about this value. 

We come to similar results in the linear approximation Sec. VI. The moments of the 
structure functions and the four-quark condensate subtract 60 MeV and 110 MeV from 
the lowest dimension value S^=270 MeV. The OPE value of the scalar self-energy m* —m 
is -140 MeV. The four-quark condensates and the moments of the structure functions add 
-140 MeV and -1-100 MeV, correspondingly. 

Turning to the role of the anomalous dimensions, we note that their inclusion into the 
moments of the structure functions lead to minor changes of several MeV of the values of 
vector and scalar self energies. Neglecting the anomalous dimensions of all the in-medium 
contributions increases the values of the vector self-energy = 230 MeV, and of the 
scalar self-energy m* — m = —140 MeV. Thus, we find for the potential energy U > if 
K = 8. However, the value of m* decreases with n, while the vector self energy practically 
does not change. We find that f/ < if k > 10, i.e. ct-^n > 55 MeV - Eq. p05|) . 

The authors of [16j carried out the detailed analysis of the nucleon self-energies 
depending on the in-medium values of the four-quark condensates. They considered 
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the QCD sum rules, based on the dispersion relations in energy go at |q| being fixed. 
The authors of P^j found that the values of the self-energies depend strongly on the 
value of scalar-scalar condensate, while the dependence on the values of the other four- 
quark expectation values appeared to be negligible small. Actually, they presented 
{M\qqqq\M) — {0\qqqq\0) = 2f{0\qq\0){N\qq\N)p, and studied the dependence of the 
nucleon parameters on the value of /. Our model calculations correspond to / = 0.14. It 
was found in [in], that the values < / < 0.3 provide the results, which are consistent 
with the nuclear phenomenology. We can deduce from Fig.l of ^H], that there values 
are m*/m = 0.65 and T^y/m = 0.28 for / = 0.14. Neglecting all the other four-quark 
condensates, we find the close values m*/m = 0.67 and H^/m = 0.25. Note, however, 
that our approach is based on the dispersion relations in another variable, i.e. in g^, with 
the relativistic pair energy s being kept fixed. (This enables us to avoid the singularities, 
connected with the excitations of medium jH]-Cni-) In our case the infiuence of the vector- 
scalar expectation value is stronger, then in ^Hl- For example, if we assume factorization 
approximation for the vector-scalar condensate, our value of the nucleon effective mass is 
about twice smaller, then the value, obtained in ^Hl- The values of the vector self-energy 
are still close in the two approaches. 

In the paper [T7] the calculations of the four-quark condensate were avoided by a 
specific choice of the function n„(g). The limits 160 MeV< <310 MeV and 0.62 GeV< 
m* <0.83 GeV have been obtained. In the work ^ the authors got rid of the four-quark 
condensates, applying the differential operators. They found the vector and scalar fields 
to be about 220 MeV and -350 MeV in the gas approximation. 

These results are consistent with each other and with the results of nuclear physics 
calculations. Various approaches in the nuclear physics studies ( see, e.g. provide 
the values between 180 MeV and 350 MeV for the vector fields, and between -200 MeV 
and -400 MeV for the scalar fields. 

There is agreement with the earlier results in some other points. The 30% reduction 
of the vector field, caused by nonlocality of the vector condensate, was found in jHl IH]. 
The strong reduction of the nucleon pole residue was obtained in jHlinilini- Also, it was 
first noted in ^] that the shift of the continuum threshold is very small. 

9 Summary 

We analyzed QCD sum rules in nuclear matter by taking into account terms of the order 
of In g^. In g^ and 1/g^ of the operator product expansion. The consistency of the lowest 
OPE terms in QCD sum rules |Hl El UHl Ull UHl HE] with the nuclear phenomenology 
was known for a long time. However the lack of information on the four-quark conden- 
sates, contributing to the terms of the order 1/g^ was the main obstacle for the further 
development of the approach. 

In this paper we studied the sum rules, treating the QCD condensates in the gas 
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approximation and included the contribution of the four-quark condensates, expressed 
through the nucleon expectation values. The latter were obtained in [20] by employing 
results of the perturbative chiral quark model \ n\ 1^ . We included also the higher 
moments of the nucleon structure functions which contribute to the terms of the order 
Ing^ and l/q^. Taking into account the four-quark condensate we included all Lorentz 
structures. 

We took into account the nonlocal structure of the vector condensate, which mani- 
fests itself through the higher moments of the structure functions. We include corrections, 
which have the smooth dependence on these moments. However, we did not include the 
nonlocal singularities in the w-channel. Such singularities correspond to the exchange in- 
teraction between the nucleon and the matter. Thus, our approach corresponds to Hartree 
description of the in-medium nucleon. The nonlocal structure of the scalar condensate 
manifests itself in the higher orders of OPE. 

Considering only the linear changes of the nucleon parameters we obtained a linear 
combination of the QCD sum rules equations in which the nucleon effective mass m* and 
the vector self-energy are the only unknown parameters. A more detailed analysis 
going beyond the linear approximation shows that this approach works well at densities 
close to the saturation value p = po- In this approach we solved the problem of expressing 
the in-medium change of the nucleon parameters in terms of the in-medium values of 
QCD condensates — Eqs. dTTIH) and (ITTT|l . 

The terms, containing the four-quark condensates provide the corrections of the order 
20 — 25% to the leading terms of the OPE of the function — Ho, which are determined 
by the local vector and scalar condensates. This is consistent with the hypothesis about 
the convergence of the OPE series. The four-quark condensates diminish the OPE value 
of the vector self -energy S^, by about 25%. The scalar self-energy m* —m is more sensitive 
to the four-quark expectation values. Inclusion of these condensates makes the OPE value 
of m* — m about 80% larger. Inclusion of the nonlocality of the vector condensate, which 
manifests itself in terms of the higher moments of the structure functions subtracts 25% 
more from the value of S^, and almost totally compensates the contribution of the four- 
quark condensates to the shift m* — m. Thus the value oi m* — m appears to be very 
close to the one, determined by the lowest orders of OPE. 

The contribution of the four-quark condensate to the vector self-energy is caused 
mainly by the vector-vector structure. The contribution to the scalar parameter m* — m 
is of more complicated origin, with the scalar-vector, scalar-scalar, vector-vector and 
axial-axial terms being numerically important. 

As it was noted earlier [^1 the QCD sum rules can be viewed as a connection 
between the exchange of uncorrelated qq pairs and the exchange of strongly correlated 
pairs with the same quantum numbers (mesons). This results in a certain connection 
between the Lorentz structures of the in-medium expectation values and of the nucleon 
propagator. In the leading orders of OPE the vector (scalar) structure of the propagator 
is determined by the vector (scalar) expectation value. The scalar-vector four-quark 
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condensate is determined mainly by the contribution which is proportional to the vector 
expectation value. On the other hand, it contributes to the scalar Lorentz structure of the 
nucleon propagator. In the meson-exchange picture such terms can be explained by the 
complicated structure of the nucleon-meson vertices. This can be instructive for model 
building of nuclear forces. 

The values of the nucleon parameters S^, and m* — m are (at least qualitatively) 
consistent with those, obtained earlier in framework of nuclear physics [32] and of QCD 
sum rules approach |7]-P, [Ilj-[I7j. The four-quark condensates, as well as the higher 
moments of the structure functions provide large contributions to the nucleon parameters. 
This future accounting for the main radiative corrections is expected to make the results 
more accurate. 

Another direction of the development of the approach is to go beyond the gas approx- 
imation. The presentation of the results, especially Eq. (Illlj) for m* enables to make the 
next step, studying the self-consistent set of equations for the nucleon effective mass and 
the quark condensates, as suggested in pUj . 
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11 Appendix A 

The contribution Xia(g) expressed by Eq. (jU^ can be obtained by direct substitution of 
Eq. into Eq. and by using the formula 

/^(ax)(6x)e^(^'^) = i 

for any vectors "a" and "6". Thus, all the contributions to the function Xia are expressed 
through the integral 

/^e^(^'^) = - ^g'Mn(-g'^) + . . . . (A2) 
Here the dots denote the terms which will be killed by the Borel transform. This leads to 

Eq. dnsi). 



(ab) + 



2{aq')(hq') 



d X 



,i{q'x) 



{Al) 
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To establish the connection between Eq. (|65p and the two terms on the rhs of Eq. 
(41) note, that the rhs of Eq. (j65j) consists of the terms of the form (see Eq. (j61|) ) 



X = J daln{Q^ + A\a))f{a) =\nQ^Jdaf{a) + J da In ^ ^ f{a). {A3) 



The first term on the rhs contains the standard logarithmic factor containing the cut, 
running to infinity. It is described by our "pole+continuum" model in a usual way. The 
second term contains a finite cut. Such terms need special treatment. The cut of the 
second term describes the singularities in the u-channel, caused by the nonlocal structure 
of the vector condensate. They correspond to the exchange terms on the rhs of the 
sum rules. We neglect such contributions, thus coming to the Hartree description of the 
nucleon in nuclear matter. 



12 Appendix B 

To obtain the coefficients of the first (factorized) terms in the rhs of Eqs. (83) and (84), 
recall that we need the expectation values of the color-antisymmetric operators 



with /i,/2 standing for the quark fiavors. The dots denote the normal ordering of the 
operators, a, a', b, b' represent the color indices. It is convenient to present 

2 1 

^aa'^bb' — ^ab'^ba' = " ^aa'^bb' ~ 7^ ^ Ka'Kb' (-^2) 

with A'' standing for the SU(3) Gell-Mann matrices Tr A^'A'^ = 25'"' . 
The factorization approximation for the quarks of different flavors is 

(M|<<(?X|M) = (M|<<|M)(M|c?X|M) (E3) 

with a, (3, 7, 5 being the Lorentz indices, and only the flrst term of the rhs of Eq. (B2) 
contributes. Using also Eq. (13) we come to Eq. (71) and (84). 

The factorization approximation formula for the quarks of the same flavor, e.g. q^^ = 
= M is 

(M|^>^'^^«riM) = (M|^>^'|M)(M|«^«r|M)- {M\uy^\M){M\uy;\M) . (54) 
Thus in the factorization approximation 

^ " {{M\uu\M)f (55) 



{M\uT^uuT^u\M) = 

16 



Tr • Tr - ^ Tr (E^E^) 
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and 

{M\Y,uT^yu-uT^XPu\M) = --Ti {T^T^){{M\uu\M)f. {B6) 
p 9 

Thus, for the factorized part of the expectation value of the color-antisymmetric operator 
rpii,uu jg = = J) 

(M|mV'mV|M) = C{{M\uu\Myf {B7) 

with 

3 V 12/ 2 V 9/ 6 ^ ^ 

Employing also Eq. (13) we come to Eq. (70) and (83). 



13 Appendix C 

Here we present for illustration the calculation of the most important contributions of the 
four-quark condensates to q structure. Using Eq. (fTSj) we find for the contribution of the 
first term of the rhs of Eq. 



= ^ ■ r [{p\i^^Up)Pp + {n\iuuUn)p,,] . (CI) 

This is equivalent to 

n' = -5M.j, (C2) 

4 

with J = J 'il}{x)il){x)d^x, while il){x) is the renormalized PCQM wave function of the 
constituent quark, normalized by the condition / ■il){x)''yo'il){x)d^x = 1. Using the value 
J = 0.54 j22|, one finds 

q^ 

The interval contribution is determined mostly by the expectation values of the oper- 
ators dr^duT^u. This happens just due to the large numerical coefficients on the rhs of 
Eq. (|7H|) . Using Eq. (ffHjl we find the contribution to be 

= (-lOa,^ - 6af - 2hl - 2bf) ^ . {C4) 

Substituting the values = -0.074^, = 0.084e[], = O.Slel = 0.06^ {eq = 
241MeV) we obtain = -0.50{el/q'^) and 

3 

n^ + n^ = 0.17 ^ . {C5) 

q2 
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A more accurate calculation, accounting for the internal contributions of the operators 
uT^uuT^u leads to the first term in the rhs of Eq. (j88p . 

The contribution to p structure is obtained in similar way. Turning to / structure, 
note that the contribution comes from the scalar- vector condensate ddu'y^u — Eq. (j85|) . 
The first ("factorized") term in the rhs provides the contribution 

= / 4^(x,r(x))e^(^^)(0|Jt/|0) p (C6) 
with e'^ix) defined by Eq. (jHD- If 0^{x) = 0^{0), we obtain 

n3 = io\dd\o) ^. 

Taking into account the dependence of 6"^ on x we actually include the higher moments 
and twists of the nucleon structure functions. Proceeding in the same way as in Sec. Ill, 
we obtain for the Borel transform of II^ 

= -Sn^mYap. (C7) 

Here 

Y = I i'-^iv) - H) - lio - miv)] m 

The first term, that is the pure local contribution, would give Y=3.0, the higher order 
contributions subtract 0.32 from this value. Thus, the factorized term would provide 
Al = 2.68. Account of the second term on the rhs of (jH^ leads to AL = 1.90. 



14 Appendix D 



Present vacuum sum rules given by Eqs. (29) and (30) in the form 
iliM^ W^) =A,+ j -^dW\ i',iM\ W^) = mAo + 



"^^0 dW 



with Ao(M^) determined by Eq. (27). In the combination £g — itlIq which is just the 
projection on the negative-parity component of the lowest state the contribution of the 
residue vanishes 

_ f (<_ ^dlV 
J \dW^ 



mil 



The condition 



d{ei - mil 



m 



dW' 



:q,I 



at = Wq means that we cannot imitate the contribution of the negative-parity pole 
of the order Aq in rhs of Eq. (D2) by changing the value of Wq. 
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15 Figure captions 



Fig.l. Solution of Eqs. (ITn^ - (ITM|l . In Fig. la the lines 1 and 2 show the Ihs of Eqs. 
f)lU2|) and ()103p for the self-energies. Line 3 shows the potential energy. The dashed lines 
show the constant values, corresponding to a„ and in rhs of these equations. The line 
in Fig. lb shows the ratio of Ihs and rhs of Eq. ()104|) . 

Fig. 2. Curves 1, 2 and 3 show the Ihs to rhs ratios of Eqs. correspondingly, at 

the values of the nucleon and continuum parameters given by Eq. ()113p . 

Fig. 3. Dependence of the solutions of Eqs. on the value of the nucleon expec- 

tation value ktv at p = po- The values of Wq, Aq are given by Eq.(33'). 

Fig. 4. Density dependence of the nucleon and continuum parameters beyond the linear 
approximation at = 8. The horizontal axis corresponds to the density of the matter, 
related to the phenomenological saturation value. 

Fig. 5. Sum-rule predictions for the dependence of the nucleon parameters m*/m and 
T^v/m on the ratio p/po at = 8. The curves correspond to the successive inclusion of 
more complicated condensates. Dashed lines: only expectation values of the operators of 
the lowest dimension g(0)7og(a; = 0) and g(0)g(a; = 0) and of the gluon operators ^6*^(0) 
are included (see Eq. lj^Uj) ). Dotted lines: local four-quark condensates are added (Eqs. 
flH3j) . (jHBj) ): sohd lines: x-dependence of the vector condensates (expressed in terms of the 
nucleon structure functions) is included. 
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